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Abstract. Spontaneous symmetry breaking in non-relativistic quantum systems 
has previously been addressed in the framework of effective field theory. Low-lying 
excitations are constructed from Nambu-Goldstone modes using symmetry arguments 
only. We extend that approach to finite systems. The approach is very general. To 
be specific, however, we consider atomic nuclei with intrinsically deformed ground 
states. The emergent symmetry breaking in such systems requires the introduction 
of additional degrees of freedom on top of the Nambu-Goldstone modes. Symmetry 
arguments suffice to construct the low-lying states of the system. In deformed nuclei 
these are vibrational modes each of which serves as band head of a rotational band. 
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1. Introduction 

In this paper we present a detailed approach towards emergent symmetry breaking 
in hnite non-relativistic quantum systems that uses concepts of effective held theory 
(EFT) m El E]. A short summary of this approach was previously presented in 
Ref. |1], but without any details. Spontaneous symmetry breaking and the use of 
EFT in inhnitely extended non-relativistic quantum systems such as ferromagnets is 
well established [HIE]. Our work is specihcally tailored for hnite systems. Examples 
are hnite Bose-Einstein condensates or BCS superconductors with a broken U{1) phase 
symmetry and molecules and nuclei that possess deformed ground states and, thus, break 
rotational invariance. In hnite systems there is no spontaneous symmetry breaking 
in the strict sense, and we speak instead of emergent symmetry breaking [7]. The 
description of such systems within an EFT requires, aside from the standard Nambu- 
Goldstone modes, additional degrees of freedom. To be specihc we address here the case 
of deformed nuclei. The additional degrees of freedom then describe rotations of the 
entire nucleus. The treatment of emergent symmetry breaking in other systems requires 
only minor modihcations, however. 

To set the stage we recall some aspects of the theory of nuclei with deformed ground 
states such as occur in rare-earth nuclei and in the actinides. Low-lying excitations 
are traditionally described phenomenologically as rotations and vibrations, both in the 
collective geometric model [8] and in the algebraic model [9]. Both models describe 
certain “leading order” aspects (or gross features) very well: Low-lying excitations are 
vibrational states that serve as band heads of rotational bands. Electromagnetic intra¬ 
band transitions are very strong, inter-band transitions are much weaker. These models 
typically fail to account quantitatively for finer details regarding “next-to-leading order” 
effects such as the change of the moment of inertia with the vibrational band head or 
the magnitude of weak E2 inter-band transitions. (Both dehciencies can be addressed 
within an effective theory mm-) 

In this paper we take the view that nuclear ground-state deformation (and, 
incidentally, pairing) are emergent phenomena. That view is supported by first ab 
initio calculations of rotational bands in light nuclei [121 [IH]. Similarly, nuclear mean- 
held calculations yield microscopic evidence for ground-state deformations na. For such 
microscopic approaches the description of emergent phenomena is a challenge, however, 
because they present multi-scale problems. Therefore, the traditional approach within 
the Wigner-Weyl (linear) realization of the underlying symmetry requires very large 
model spaces to describe the emergent symmetry breaking. 

The approach using an EFT is based entirely upon symmetry arguments. The 
emergent broken symmetry is treated using the nonlinear Nambu-Goldstone realization 
of the symmetry (as opposed to the linear Wigner-Weyl realization) plus those additional 
degrees of freedom that account for rotations. The combined treatment of these degrees 
of freedom is the novel technical aspect of our work. The approach has the advantage 
of being model independent. A controlled expansion in terms of well-dehned small 
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parameters generates higher-order terms is a systematic fashion. These may serve as 
useful guides in the phenomenological models. 

Aside from presenting details of the EFT approach outlined in Ref. |1], it is the 
purpose of this paper to show how such higher-order terms are generated and how they 
influence results of lower order. Moreover we establish the connection of our EFT with 
the technically simpler effective theory of deformed nuclei developed in Ref. [15]. Our 
results show that the effective theory is based upon a solid field-theoretical foundation. 

2. General Approach 

We consider a deformed liquid drop with a space-fixed center of mass and with axial 
symmetry about the body-fixed 2 ;'-axis. In the present Section we introduce the Nambu- 
Goldstone modes and in addition the time-dependent modes that describe rotational 
motion. We warn the reader ahead of time that in Section [3] we switch to another 
parameterization of these modes. The parameterization used in the present Section 
is physically transparent and easy to justify. It has the drawback that it leads to 
analytically cumbersome expressions. The parameterization used in Section [3] is easy to 
handle but based on a different physical picture. In order to exhibit the general scheme 
we follow the approach of the present Section up to the construction of the classical field 
theory and its symmetries. The developments of Section [3] then run in parallel to these 
developments. 

2.1. Nambu-Goldstone Modes 

We denote the Cartesian coordinates in the space-fixed system S (the body-fixed system 
S') by {x, y, z} (by {x', y', z'j, respectively). The nuclear ground state is invariant under 
50(2) rotations about the body-fixed 5-axis while 50(3) symmetry is broken by the 
deformation. The Nambu-Goldstone modes parameterize the two-dimensional coset 
space 50(3)/50(2) and depend on two complex fields 7r±. In the body-fixed system, a 
mass element of the liquid drop has spherical coordinates (r, 6*,0). The classical fields 
TT± depend on these dynamical variables and on time t. We neglect the r-dependence of 
7r± because radial vibrations of the liquid drop are expected to have higher excitation 
energies than surface vibrations. We denote the remaining coordinates 6, (p and time t 
jointly by with p = 1,2, 3. The fields 7i± generate the local rotation 

U{p) = exp{-in-{p)P+ - z7r+(p)P_} . (1) 

We follow Ref. [16] and denote the total angular momentum operator in the body-fixed 
system by P and in the space-fixed system by J. Then P^', Pp, Pz' are the components 
of P in the body-fixed system, and 

P± = z-—j^(^P^i ± iPy') . (2) 

Although we use the expression “angular momentum operator” for P and J, the three 
components of these operators represent in the present context only the three generators 
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of rotations in the body-fixed and in the space-fixed system, respectively. With 7r_|_ = tt^, 
the transformation U is nnitary. For a physical interpretation of the helds 7r± we write 

7r± = ^a;exp{±<} (3) 

with uj and ( as the new real dynamical fields so that Eq. ([T]) becomes 

U = exp{—za;[cos CF’x'+ sin ^Py]} . (4) 

We interpret Eq. (|1]) geometrically. In the body-fixed system, rotations are aronnd an 
axis perpendicular to the symmetry axis. The azimuthal angle of the rotation axis is 
denoted as The angle u is the angle of rotation about that axis. Replacing ^ by C -|- vr 
we change the orientation of the axis of rotation or, equivalently, the direction of rotation 
(clockwise rotation —)■ counter-clockwise rotation). Clockwise rotation by the angle uj 
is equivalent to counter-clockwise rotation by the angle 27 r — uj. Therefore, we might 
conhne the angle cu to 0 < a; < tt. For reasons that will become apparent later we do not 
adopt that choice here so that the ranges of the angles ( and u are 0 < cu, C < 27r. We 
expect the effective Lagrangian to depend on trigonometric functions of u and ( only. 
The fields uj and ( are functions of the dynamical variables jji = 1,2,3. Inspection 
shows that a pure time dependence of u and ( describes the rotation of the deformed 
nucleus as a whole while a genuine dependence of u and ( on the angles 9 and (f accounts 
for surface vibrations, i.e., local dislocations of the constituents of the deformed drop 
(fluid elements or nucleons, as the case may be). Therefore, we write 

u = uo{t) +Ui{9,(j),t) , 

C =Co{t) + Ci{0A,t) , (5) 

with the understanding that Ui and Ci possess a non-trivial dependence on at least 
one of the angles 9, 0 (so that not all partial derivatives with respect to these angles 
vanish identically). We expect (and verify later) that the functions uji and Ci define the 
true Nambu-Goldstone modes of the deformed nucleus that describe surface vibrations. 
For these we will use a small-amplitude approximation. The purely time-dependent 
functions uJo{t) and Co(^)) on the other hand, are alien to the standard approach to 
spontaneously broken symmetry. They represent the novel element that accounts for 
the finite size of the system and describe rotational motion. For these functions we 
cannot use the small-amplitude approximation. This situation is technically similar to 
hnite-volume EFTs of quantum chromodynamics [ni[l8] and magnets [19] . Here - and 
in contrast to the genuine finite systems we consider ~ the finite volume stems from 
limitations of computational resources and obscures the physics of the infinite system 
one is interested in. In both cases, the purely time-dpendent “zero mode” exhibits large 
fluctuations and has to be treated separately. 

We follow Ref. [5] and define the Nambu-Goldstone modes and of our problem 
by writing 


U*(p)td^U(p) = 0-P+ + o+P_ + . 


( 6 ) 
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Here = d/dp^j, for p = 1,2, 3. The coefficients and may alternatively be written 
as functions of the helds or of the helds u, (. The effective Lagrangian is obtained 
by forming combinations of these coefficients that are invariant both with respect to the 
group operations and with respect to space rotations. 

We express the Nambu-Goldstone modes in terms of u and (. The calculation is 
lengthy but straightforward. We use the left-hand side of Eq. ([6]), the form (jl]) for U, 
and a Taylor expansion for U. We hnd 

% = (l/2)exp{-<}[a/,a;-z(a^C) sinca] , 

= (l/2)exp{+iC}[5MW-l-i(a/,C) sinw] , 

<'= (^/^C)(l-cosa;) . (7) 

2.2. Invariants 

To determine the behavior of and under group operations we consider the action 
of a hxed element g (independent of the of the coset space on U. Equation ([T]) 
implies that under the action ol g, U changes nonlinearly, 

U^[gU]h^{g,U) . (8) 

The product [gU] lies within the coset space SO{3)/SO{2) while 

h{g,U) = exp{z^{g,U)P,^} (9) 

is an element of 50(2). The function T depends on the parameters characterizing both 
g and the transformation U. Thus, 

U^{p)td,U{p) ^ hU^{p)t[d,U{p)]h^ + d,^{g, U)P,, (10) 

and 

a^P+ -I- a'^P- + Pz' h{a^P+ + a'jfP- + Pz')h^ + U)Pz' .(11) 
We recall the commutation relations 

[4, Jy] = iJz (cyclic) , 

[Px',Py'] = -iPz' (cyclic) , 

[Jk,Pi'] = 0 (aW k,l = x,y,z) . (12) 

which imply 

[Jz^ J+\ J+ ) ['Jzi J—] J— 5 [<7+, j—] Jz 5 

[Pz^, P+] = -P+ , [Pz>, P-] = P- , [P+, P-] = Pz' (13) 

and PfIP± = P±{Pz' T 1)”^- To calculate the right-hand side of expression flTT]) we use 
these for P and the expansion of h in a Taylor series in Pz'. We obtain 

% exp{-iT}a; , 
a+ ^ exp{iT}a+ , 

+ dy^ig, U) . 


( 14 ) 
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Therefore, the terms of lowest order in the a^, that are invariant under group 
operations are 0 


a+Oj, and . (15) 

Here /i, z/ = 1, 2, 3. The expressions (IT^ are obtained by treating the derivatives both 
with respect to the angles 6, 0 and with respect to time as though they were derivatives 
with respect to external variables. While that is appropriate for time t (/r = 3), it is 
not for the angles 0 (/i = 1, 2) because these change under rotations. The effective 
Lagrangian must be invariant with respect to such rotations. In actually constructing 
the invariants we demand axial symmetry (see Section [3]) and are guided by the analogy 
to coordinate transformations in Euclidean space in three dimensions. Here we would 
interpret dk with k = x,y,z as one component of the vector V, use the invariance of 
the scalar product of two vectors, and obtain the invariants Va'^'Va~. On the unit 
sphere we analogously replace V by the vector L of orbital angular momentum with 
components 


r ■■ 1 ^ ■ 1 n ^ 

Lx' = + «sm 0— + i cos 0 cot a-— 
oO od) 


d 


d 


Lyi = — i cos 0— + i sin 0 cot 9— , 


09 


00 


r ■ ^ 


The lowest-order axially symmetric invariants formed from are then 


( 16 ) 


d-jx'a J-jx'a I d-jy^a djyia , 

Lz'a'^Lzia~ , 

dta'^dta~ . (17) 


Expressions like Lx'a~^ are here understood as linear combinations of terms (read as 
partial derivatives with respect to py and explicitly given in Eqs. (I7j)) with coefficients 
given in the hrst of Eqs. flT6|l . For simplicity, we employ in the present Section only the 
spherical invariant 

Lx'O^Lx'O T Lyio^Lyia T L^'o'^L^'O , (18) 

and use the more general axially symmetric invariants in Section |3l 

The substitution rules flTdll apply in both the body-hxed and the space-hxed 
coordinate system, albeit with different dehnitions of the function T. Using this fact 
and the well-known behavior of the components of L under rotations we easily conhrm 
that the form (ITSll is indeed rotationally invariant. 

The eigenstates of the nuclear Hamiltonian are (almost exactly) eigenstates of the 
parity operator. But that operator cannot be written in terms of rotations. Indeed, 
in three dimensions the determinant of the matrix describing parity inversion equals 
minus one while the determinants of all rotation matrices are equal to plus one because 
these connect continuously to the unit matrix. Therefore, we cannot incorporate parity 
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conservation into the present approach. We retnrn to that point in Section 14.21 The 
situation for time reversal differs. To be invariant under time reversal the effective 
Lagrangian must contain only even derivatives with respect to time. That constraint 
has been taken into account in expressions (fT7|) . The constraint is appropriate for even- 
even nuclei with ground-state spin zero. Another discrete symmetry (the 7^-parity of 
Ref. |8]) plays a role for the quantized version of the theory and is treated in Section 
too. 

The explicit form of the invariants is obtained by using Eqs. ([7]) in expressions ca 
and flT8|) . 

Cl = 2afaf = ^(d;‘^ + C^sin^o;) , 

3 .. 

£2 = 2^(Lfca+)(Lfca") =-(^(Lcn)^ + sin^o; (LC)^) • (19) 

k=l 

The dot indicates the time derivative. 


2.3. Classical Field Theory 

The classical held theory is obtained by writing the effective Lagrangian density £ as a 
linear combination of the two invariants (IT^ . 

C = CCi + DC 2 . (20) 


The coefficients C and D are the parameters of the theory and are determined by a ht 
to the data. The total Lagrangian L is obtained by integrating C over the dynamical 
angles 0,0, 


r-27r 


L= dE C = 


dvr 


d6 sin 6 


d0 C . 


( 21 ) 


The normalization is chosen such that for C = 1 the integral in Eq. fl2T|) yields unity. 
The factor sin 6 in the measure fITTD takes account of the fact that 6, 0 dehne points 
on the surface of a sphere so that we deal with curvilinear coordinates. We show in 
Appendix 1 that in such coordinates the equations of motion are 




1^9, 


9(sin 0£)\ 

d{d^u}) ) 

d{sm 6C) 

9{d,C) 


c?(sin 9C) 
doj 

9(sin 6C) 

dC 


( 22 ) 


We recall that /r = 1, 2 stands for the angles 0, 0 and /r = 3 for the time t, and that 
both oj and ( are functions of 0, 0, t. Eqs. fl22|) . fl20|) . and flT^ constitute the non-linear 
equations of motion for the two classical helds u and (. 

In constructing C we have imposed rotational invariance. Therefore, we expect that 
the total angular momentum of the system is conserved. Using the Noether theorem as 
in Ref. [15], we now construct three constants of the motion. These correspond to the 
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three components of angular momentum (not to be confused with the three components 
of the rotation operator in Eqs. flT^ i. With k = x',y',z' we consider the rotation 


r = ex 


SXtPi} 


(23) 


by inhnitesimally small angles Sxk about the three body-fixed axes. The ensuing 
inhnitesimal changes of u and of ( are denoted by 6u and 6(, those of C by 6C. A 
straightforward calculation yields 




d ( dC 


k=l 


^ 9/r \d{d^u 


-Alik + 


dC 


dm) 


M 


2k 


(24) 


With 5uj = 5qi and 5(f = 5q2 the matrix M is dehned by 5qi, = M^k^Xk and 
explicitly given in Appendix 2. Rotational invariance implies = 0 for every choice of 
Sxk with k = x', y', z'. That implies the vanishing of each of the three Noether currents 
k = x',y'^ z' in big curly brackets in Eq. fl2Tl) . 

To obtain the constants of the motion we consider an inhnitesimal change of the 
total Lagrangian given by 

6L = [ ^^6(smeC) , (25) 

J dvr 


with SC dehned in Eq. (l24|l . The ehective Lagrangian density C is periodic in the angle 
0. With respect to 6, the same statement holds for /Isin^. Partial integration then 
shows that the derivative terms with respect to 6 and 0 in Eq. fl2Tll do not contribute 
to 5L so that 5L is given by a pure time derivative. The vanishing of 5L for any choice 
of the angles 6xk then implies the existence of three constants of the motion. These are 
the three components Qk of the total angular momentum of the system. With M given 
in Eq. (11181) of Appendix 2 and the integration measure defined in Eq. fl2TD these are 


Qx' 


Qy' 


Qz' 



( 26 ) 


2.4- Power Counting 

Without much justihcation given, we have constructed the effective Lagrangian in 
Eqs. fl20|) and flT^ from the lowest-order invariants in the Nambu-Goldstone modes. 
Further progress hinges on the identification of the energy scales that govern our problem 
and of the associated powers of higher-order terms in the Nambu-Goldstone modes. 
The relevant scales are the energy scale of rotational motion, the energy scale G of 
vibrational motion, and the cutoff parameter A beyond which other modes like fermionic 
excitations play a role. In our approach it is assumed that ^ -C G <C A. 
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To identify the various contributions to the effective Lagrangian in Eq. 0201) and to 
the equations of motion fl22|) . we use the decomposition ([5]) and an expansion of Ui and 
of Cl in terms of spherical harmonics, 

CO 

UJl = EE 

L=2 

OD 

c> = EE CL,{t)YL,{e,<f) . (27) 

L=2 fi 

The term with L = 1 describes center-of-mass motion and is suppressed. As mentioned 
earlier we take rotational motion fully into account by treating uo and Co without any 
approximation. We use a small-amplitude approximation for the surface vibrations 
described by ui and Ci- A naive hrst approach would then consist in expanding C in 
powers of Ui and of Ci and in keeping only terms up to second order (but terms of all 
orders in uq and in Co)- That would be in line with the treatment of ferromagnets and 
paramagnets in Ref. [5] where in the leading-order effective Lagrangian only terms up 
to second order in the helds are kept. However, because of the presence of rotational 
motion that procedure does not fully apply in our case. To see that we must analyze 
the contributions to C in some detail. We use the symbol ~ to dehne the relevant order, 
and we drop the magnetic quantum numbers on ul and Cl- 

The ratios Cjq/ujq ~ C ci-nd Co/Co ~ C ci-re governed by the energy scale of rotational 
motion. The range of the variables ujq and Co is of order unity. Therefore, (ho ~ C ^md 
(^0 ~ C- Inserting Eqs. ([5]) into Eqs. fl2^ we obtain the terms (C'/2)(ho and {C/2)Q. 
We show below that these describe rotational motion. Therefore, {C/2)6jq ~ C ^md 
{C/2)Q ~ c- Together with (ho C and Co r'U C that implies C ~ 1/C which is consistent 
with the interpretation of C as moment of inertia. The ratios uJi/ujl ~ H and Cl/Cl ~ 
with L >2 are similarly governed by the energy scale H of vibrational motion. Inserting 
Eqs. (Ej) into Eqs. fl2^ we obtain for L > 2 the terms (C/2)(h£ and {C/2)Cfl. These 
describe vibrational motion. Therefore (C/2)(h| ~ H and [C/2)Cff ~ H. Together with 
C ~ 1/C these relations imply col,Cl ~ \/^ and, together with col/ojl Cl/Cl 
also UJl Cl \/C/kl for L > 2. These relations are used below when we 

expand C in powers oi ul and Cl- The terms in £2 (see Eqs. flT^ ) describe vibrational 
motion so we have {D/2) {Lul)^ ~ fl and {D/2){L(l)‘^ ~ fl. 

In our approach the operator L is dimensionless (see Eqs. ffT6l) L Formally, however, 
L plays the same role as the momentum in the theory of Ref. [5] where only small 
momenta are kept for small energies. Physically we analogously expect that only 
the small eigenvalues of the operator are relevant for the low-energy part of the 
spectrum in our case. Therefore we formally attach to L the scale H. Together with 
Cl ^JC/Et that gives D ~ 1/C. We summarize these assumptions and results by 
writing 


c < , 

1^0, Co ~ 1 , 
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ujl, Cl ~ Vc/^ -C 1 for L > 2 , 
'^o,Co ~ C , 


dlo, Co 




(Al, Cl ~ 
Al, Cl ~ 


C 

D 

L 


for L > 2 , 
for L > 2 , 

1 

C ’ 

1 

c ’ 

n. 


(28) 


These considerations imply the following rnle for an approximate treatment of the 
problem. The effective Lagrangian C in Eqs. fl20l [T^ contains terms of order hi that 
describe nnclear surface vibrations and terms of order C that describe rotational motion. 
In expanding C in powers of ul and Cl we must, therefore, keep terms of orders hi, yTlC 
and C- We omit terms of order \JC/Et or less. The resulting approximate expressions 
are 


C ( 

CCi ^ — (^(cuo + Aif + Co sin^ uq 


2CoCi 


sin^ uq + oji sin 2u;q 


+ Ci 


sin^ uo + ui sin 2a;o + cos 2a;o 


DCo 


D 


(Luf + {LCf 


sin^ uq + oji sin 2u)o 


-f- cos 2 cUq 


(29) 


These developments show which terms to keep in the effective Lagrangian. 


2.5. Purely Rotational Motion 

To ascertain consistency of our arguments we consider the case of nuclear rotation 
without surface vibrations. We accordingly assume that uj and C depend only upon 
time, so that in Eqs. ([5]) we have cui = 0 = Ci- For purely time-dependent fields 
classical field theory changes into classical mechanics. The treatment becomes very 
similar to that of Ref. [12]. The effective Lagrangian is 

^ = ^{aI + Co sin^ cjq) , 

and the canonical momenta are 
dC 


( 30 ) 
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The Hamiltonian H is given by 


H = 


2C 


+ 


1 2 
2 , , ’^C 


Sin ujq 


The three components Qk of angular momentum are given by Eqs. 

Qk = TTujMik + 'n-(^M2k ■ 

We use Eq. flll8|) of Appendix 2 for M and hnd 


'^Ql = T^l + - 

( ^ Cl 


1 


2 

sm ojq ^ 


(32) 

which now read 

(33) 

(34) 


That shows that the Hamiltonian 
momentum, 


1 




2 

k ■ 


is proportional to the square of the total angular 


(35) 


In other words, we obtain the classical theory of the rotating top. The constant C is 
the moment of inertia. 


3. Another Parameterization 

We have actually carried the approach of Section [2] further, deriving the Hamiltonian 
and quantizing it. The resulting equations are difficult to interpret, however. They 
do not display in an obvious fashion what is expected on physical grounds: Harmonic 
vibrational motion of the variables ul and C,l- As shown in Appendix 3, these difficulties 
have to do with the non-Cartesian form of the measure dE in Eq. fl?I|) . That is why we 
now introduce another parameterization of the matrix U defined in Eq. ([T]). We proceed 
in close analogy to Section [2l 


3.1. Nambu-Goldstone Modes 


We use the space-fixed system, and we parameterize the matrix U in product form. 


U = g{C,u}) u{x,y) , 

g{C,uj)= exp |-iC(t)^ 2 }exp |-iu;(t)4| , 

u{x,y) = exp ^^-ixJ^ - iyjy^ . (36) 

The purely time-dependent variables u and ( describe rotations of the finite system, 
similarly to the variables uqXo introduced in Eqs. ([5]). As in Section [2] we choose the 
ranges as 0 < ca, (^ < 27r. This is convenient for Section 14.11 With 6 and (j) as defined 
in Section [2|, the fields x = x{9, 0, t) and y = y{9, 0, t) play the role of the fields Ui and 
Cl defined in Eqs. ([5]). They describe the small-amplitude vibrations of the liquid drop. 
To exclude the possibility that x and y induce a global rotation of the entire drop we 
request 


J dE x{9, (/), t) = 0 


dE y{9,((),t) ■ 


(37) 
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We use the definition that U = gu acts onto objects to the right. Thus, the local 
vibrations induced by the field u are followed by a global rotation g of the entire drop. 

We show in Appendix 3 that the parameterization of U in terms of the variables uj 
and C in Eqs. (jl]) and ([5]) and the one introduced in Eq. fl36|) are completely equivalent. 
Why then did we not start from the outset with the parameterization fl36ll ? As shown 
in Section [21 the parameterization used in Eqs. (jll) and ([5]) can be justified physically in 
a convincing manner. Moreover, it is tailored after the standard approach to symmetry 
breaking in non-relativistic systems. The advantage of the new parameterization is 
that it treats the rotational degrees of freedom separately while the parameterization 
in Eqs. (jl]) and (JSl) treats the rotational mode and the vibrational modes on an equal 
footing. It has an alternative physical interpretation. When acting from right to left, u 
induces a small-amplitude dislocation of a nucleon (or volume element) at (0, 0) in the 
axially-symmetric nucleus (whose symmetry axis is the z axis), while g then rotates the 
entire nucleus. 

Power counting as in Section 12.41 shows that 

a;,C ~0(1), 
kl, \y\ ~ < 1 , 

x,y ~ . (38) 

The parameter e helps to identify (and omit) higher powers of x and ?/, consistent with 
a focus on small-amplitude harmonic surface vibrations. A physical interpretation of 
this parameter is given at the end of Section 13.31 below. From here on the development 
is similar to that of Section [2l 

In analogy to Eq. ([6]) we define 

U-Hd,U = afJ, + alJy + afJ, . (39) 

As in Section [2] the symbol with /i = 1,2,3 stands for the partial derivatives with 
respect to the angles 0, 0 and time t while in the present Section d„ with z/ = 1, 2 stands 
for the derivatives with respect to the angles only. To work out the Nambu-Goldstone 
modes explicitly we use 

g~^dtg = — i (^—(smuJx + ujJy + ( cosuJz') , 

u~^dtu = -i(^x + ^{xy- yx)^ 

-i {y - - yx)) Jy 

-i (^{yx-xy)^ Jz , 

and, from the Baker-Campbell-Haussdorf expansion, 

U-^d^U = u-^dyu , 

U~^dtU = u~^dtu + u~^ {g~^dtg) u 


( 40 ) 
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= — i [x — ( sin N ^ {xy — yx) — y( cos ui ) Jx 


X 


— i yy + ui — —{xy — yx) + xC coscuj Jy 

— i{ - {yx — xy) + C cos a; — xcj — ?/C sin a; ) Jz + .. 


(41) 


Here and in what follows, the dots indicate terms of higher order in e. From Eqs. 
and fl4T]) we obtain 

y 

af = X -\— {xy — yx) — ( sinoj — y( cos a; + ... , 

6 

of = y - {xy — yx) + ci; + coscn + ... , 

6 

af = — - {xy — yx) + ( cos oj — y( sin oj — xoj + ... , 

and 


(42) 


al = d^x + ^{xd„y - yd„x) + ... , 

D 

< = duV - -^{^dyV - ydyx) + . . . , 

0 

< = - ]^{xd^y - ydyx) + ... . (43) 

Eqs. fl42|) and fl43l) give the lowest-order contributions to the Nambu-Goldstone modes 
for the parameterization fl36l) . 


3.2. Invariants 

As in Section 12.21 we build the effective Lagrangian upon invariants constructed from 
the Nambu-Goldstone modes. To this end we need to determine the behavior of these 
modes under transformations. We consider a rotation r about infinitesimal angles 6xk 
around the space-hxed k = x,y,z axes. We use Eq. fl36|l for U. With 

rg{C,^) = 9{C\^')h{'y') , Hi) = ewHiH} , (44) 

and C'yX)', 7 ' given below we have 

rU = rg{C,uj)u{x,y) = g{C,uj')h{i)u = g{C,uj') [h{-f')uH{H)] h{i) = U'h .(45) 

The last of equations fj45ll defines U' = U{(',u',x',y') as an element of the coset space 
S'0(3)/S'0(2). We accordingly write 

U'-HdyU' = {aD'Jx + {aD'Jy + {aD'Jz . (46) 

Proceeding as in Section 12.21 we have 

U-^dyU ^ {rU)-^dy{rU) 

= {U'h)-^dy{u'h) 

= h-\U')-\dyU')h+h-\dyh) . 


( 47 ) 
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From Eqs. (147114^ and Eq. 


and 


we obtain 



cos 7 
sin 7 ' 



cos 7 


(<)' = < + 


(48) 


(49) 


We have used that 7 ' as given in Eq. fl52l) below is a function of t only. 

It remains to work out the relation between the variables ca', x', y' and (f, u, x, y. 
A calculation similar to that of Appendix 2 shows that g{(,uj) transforms into g{(',uj') 
with = C + and co' = oj + 6u where 

^ ^Xx ^ 

6xy ■ (50) 

V 

According to Eq. fl45|l the matrix u transforms under the action of r into h{'j')uh^{'j'). 
That transformation differs from that of Eq. (IHD because of the prefactor g in the 
definition of U in Eqs. (1361) . Therefore, x and y transform into x' and y' according 
to 



cot u cos 
-sinC 


cot ca sin 1 

cos 0 


X 


cos 7 ' — sin 7 ' 
sin 7 ' cos 7 ' 



where 


, cos C e- sin C e- 
7 = — 


(51) 


(52) 


sm uj sm uj 

Eq. flSTl) shows that is invariant under rotations. Moreover, since 7 ' depends 

on time, under rotations the four quantities x, y, x, ij are transformed into linear 
combinations of x', ?/', x', y'. 

We are now ready to construct the invariants. We begin with the time derivatives 
in Eqs. (1421) . Eq. (1481) shows that (af)^ + (a^is invariant. We use the power counting 
of Eqs. (138|) (see also Section 13.41 below) and drop terms of order ffe and with 
k > 3/2. We also omit terms linear in x, y, x, or ij as these vanish upon integration over 
6 and 0, see Eqs. (137)) . That gives 

[afY + {afY An^ oj-\-x^+y‘^+2{xy—yx)C,cosoj — ^{xy—yxY .(53) 

o 


The invariant form (15^ is the sum of three homogeneous polynomials of orders zero, 
two, and four, respectively, in the variables x, ?/, and their time derivatives. Under 
rotations, each of these is transformed into another homogeneous polynomial of the 
same order, see the text below Eq. (I52l) ). Invariance of the form (l53l) implies that each 
of the said polynomials is invariant by itself. Hence the invariants are 

E\a — ^ + C ^ ) 

= x^ + ?/^ + 2 (x?/ — yx)C, cos a; , 

■Cic = m - yxf . 


(54) 
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The additional invariant 

+ 2/^) [x‘^ + y‘^ + 2(a;y — yx)( cos u] (55) 

is obtained by multiplying Cib with the invariant + y‘^). The invariant Cid is of the 
same order as Cic- 

We turn to the invariants constructed from the derivatives with respect to the angles 
6^,0 in Eqs. (03]). We conhne ourselves to terms of up to fourth order in x and y and 
their derivatives. Eq. f08|) shows that for all u = 6,(j) the form + (a)))^ is invariant, 
and so are al and alaf. Forming suitable linear combinations of these and multiplying 
with the additional invariant (x^ + y'^) we hnd the invariants 

C 2 a = {Lxf + {Lyf , 

C'2a' = {LzXY + {L^yY ) 

C 2 b = {xLy - yLxf , 

Ac ={x^ + y‘^){{Lxf + {Lyf) . (56) 

The construction of the invariants in Eqs. (I5T|) to fl56|) is based on the fact that under 

rotations, u transforms into huhf This feature does not apply in the absence of any 
rotational motion, i.e., for = 1 or cu = 0 = C- If is easy to see that in that case 

we would have rU = ru = u{x',y')h where x',y' are nonlinear functions of x,y. The 

argument shows why the invariants in Eqs. fl5T|) to fl5^ occur specihcally in the case of 
rotational motion but differ in cases like ferromagnetism or paramagnetism where all 
modes considered are true Nambu-Goldstone modes. It also shows that care is needed 
when considering the limit of an inhnitely large moment of inertia: In the framework of 
the present formalism, that limit differs from the one where rotational motion is ruled 
out from the outset. 

In the construction of the invariants, time-reversal invariance has been taken into 
account in the same manner as in Section 12.21 When we apply our formalism to atomic 
nuclei, an additional symmetry (the 7^-symmetry) comes into play. That symmetry 
matters for the quantized version of the theory and is, therefore, deferred to Section 14.21 

3.3. Classical Field Theory 

As in Section 12.31 the effective Lagrangian L is given in terms of an arbitrary linear 
combination of the invariants constructed in Section 13.21 and involves eight constants 
Gj, i = a,b,c,d and Di, i = a,a',b,c that must be determined by a £t to data. In 
obvious notation we have 

L = Li + L2 = J C 

i=a,b,c,d i=a,a',b,c 

The integration over angles is dehned in Eq. (ETj). 


( 57 ) 
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In slight difference to Eqs. fl27p we expand the real variable x in two ways, either 
in spherical harmonics = (—or in terms of the real orthonormal functions 

f 71 ’ /i > 0 , 

/i = 0 , (58) 

[ 7^ ) /i < 0 . 

We note that the functions do not form the components of a spherical tensor. We 
write 

CO L 

X = ^ ^ XL^,ZL^, 

L=‘l fi=—L 

CO L 

= 5^ , (59) 

L=2^=—L 

and correspondingly for the real variables y, x, y. As in Eqs. fl27)) . contributions with 
L = 0 and L = 1 are excluded. The coefficients are real. For the complex coefficients 
xl^ we have xf^ = (—For every value of L the coefficients x^^ and are 
linearly related in an obvious way. The hrst of Eqs. (I5^ is useful for quantization. The 
second is more useful when the calculation requires angular-momentum algebra. 

Using the hrst Eq. fl59|) and carrying out the integration over angles, we obtain for 
the total Lagrangian in Eq. fl57)l 

L = ^ (72 + C^sin^cn) + y {x% + 

L n 

+ CbC cos a; EE {xLEijLfi yLfiXL^i) 

L fj, 

~ 2 {Y>aL{L + 1) + Daiy^) {x\^ + yf^) . (60) 

Lfi 


We have restricted ourselves to terms up to and including the orders 0{Et), 0{eVL), and 
0{ff). The canonical momenta are 


Pu 



dC „ dC y dC 


(61) 


In analogy to the hrst Eq. fl59l) we dehne the real function 

CO L 

= < ( 62 ) 

L =2 ti=—L 


and correspondingly for Py{6, cf). 

We use the Noether theorem as in Section 12.31 and hnd the conserved quantities 
{k = x,y,z) 


Qk= dE 


6u 6( 6x 6y 

+ ^Pc + ^Px + ^Py 

oXk oXk oXk oXk 


(63) 










Effective Field Theory of Emergent Symmetry Breaking in Deformed Atomic Nuclei 17 


We identify the Qk with the three components of angular momentum. Explicitly we 
obtain from Eq. flSOjl and from the differential form of Eq. flSTj) 


Qx 

Qy 

Q. 


. ^ ^ COsC 

— sm C — Pc cot oj cos C H- K 

sino; 

sin 

p^ cos C ~ Pc cot uj sin C H- K , 

sino; 

Pc ■ 


(64) 


Here 

K = J dE {xpy - yp^) . (65) 

The terms in Eqs. (1M|1 that do not involve the factor K correspond to the angular 
momentum of the rigid rotor. That is shown below and is analogous to Section [231 The 
integral K over solid angle in Eq. fl65l) is the angular momentum of the two-dimensional 
oscillators that describe the surface vibrations. That is easily shown by applying the 
Noether theorem to / dE {x^ + ?/^). As remarked below Eq. (l52|l . that expression is 
invariant under SO{2) transformations of {x,y). The conserved quantity associated 
with this invariance is K. The square of the total angular momentum is 

= pI+ . ^2 (Pc ~ ^KPC cosca + 77^) . (66) 

sm cj 


3.4- Power Counting 


Given the full Lagrangian in Eq. fl60|) we can now complete the arguments leading to 
the relations (l38ll . To identify the terms that are kept we use arguments similar to the 
ones in Section 12.41 In addition to Eqs. (138|) we assume 


Da, Da' ~ G/e , 

C'a ~ r', 

which implies 

CaEla ~ ^ , 

CbEib ~ D , 

CcE\c, CdEid eSI , 

and 


(67) 


( 68 ) 


Puj,PC ~ 1 , 

Px,Py ^ ■ (69) 

If we were to scale Cc ~ (e^G)”^ we would hnd CcEic ~ G, and that would be unexpected 
(or “unnatural”) for a term with such a high power in the coordinates x, y. We cannot 
completely rule out this possibility, however, as effective held theories with unnaturally 
large scale do exist. As an example we mention the pion-less nuclear effective held theory 
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for systems with a large scattering length [201 E]- The large values of Px,Py implied by 
the relations fl69l) do not seem natural but cannot be avoided if we insist on x,y ^ 

To illuminate the role of the dimensionless parameter e we consider the limit of 
an inhnite system (with infinite moment of inertia Ca, or ^ 0). The fields C, and 

u become static and have constant values that depend on the orientation of the rotor. 
Spontaneous symmetry breaking gives rise to the Nambu-Goldstone helds x and y that 
describe the low-energy modes at the scale G. The effective field theory breaks down 
at the scale A S> G. The assumption e: ~ G/A <C 1 takes account of this fact and is 
similarly used in other non-relativistic applications of effective field theory mu- Terms 
of higher order in the time derivatives of a; or y that were suppressed in Eq. fl5^ would be 
of order with k > 3/2 and, thus, of higher order in G/A. A similar power counting 
for the spatial derivatives results from the replacement L —)■ QL as in Section 12.41 

In the opposite case where A —)■ cxo but where ^ <C G is finite (i.e., differs from 
zero), the terms of order oc 1/A disappear. The Hamiltonian describes a rotor 
coupled to a set of oscillators. A problem occurs once the excitation energy is so large 
that the amplitudes x, y are of order unity and compete with the finite rotations of the 
top. A distinction between the two types of modes is then no longer meaningful, and 
spontaneous symmetry breaking does not give an adequate description of the system. 


3.5. Effective Hamiltonian 


The effective Hamiltonian H is obtained from the effective Lagrangian L via a Legendre 
transformation. To perform that transformation we write the kinetic part Li of L 
in a form that displays its bilinear dependence on the velocities. We define the 
inhnite-dimensional velocity vector V'^ = {(/,ci;,where L = 2,3,... and 
p = L,L — 1,...,—L and write Li = {1/2)V'^GV. The matrix G is easily found 
from Eq. fl60l) . We analogously define the vector of momenta = {pc,,Pu}iPfy,pffj}- 
Then the effective classical Hamiltonian is 


H = \p^G-^P + i {DaL{L + 1) + Da'p^) 


X 


LfL 


+ y\fi) 


(70) 


Lfi 


To calculate the inverse G~^ we write G = Gq + Gi, where Gq is the diagonal part of G, 
and use perturbation theory in Gi so that G~^ = Gg ^ —Gg ^GiGg ^ + Go ^GiGq ^GiGq 
.... Keeping only terms up to and including the orders G(H), 0{ekl)^ and G(^) we obtain 


H = 


1 


^ {{plfjf + + {DaL{L + 1 ) + Da>p^) {x\^ + y\^ 


LfL 

1 


pI + 


sin^o; 


+ 2Kp(^ cos u + cos^ uf) 


2 ^ {{PliX + {PhyY) + {DaL{L + 1 ) + Da'P^) {x\^ + yf^) 

Lfi *- 


( 71 ) 
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We have used Eq. (IM]) and [see Eq. fl^ ] 

^ = ■ (72) 

Lfi 

In the Hamiltonian fTfll) . the Nambu-Goldstone modes undergo harmonic vibrations. 
These are coupled via to a rigid rotor. The vibrations are of order 0{VL), while the 
rotations are of order 0{f). 

We mention in passing that an alternative (and simpler) form of the Legendre 
transformation seems to exist. Instead of defining the momenta as in Eqs. fl6T|) and using 
Eq. (1701) . one might dehne the canonical momenta px{0,(j)) and Py{0,(j)) as functional 
derivatives of Jf^CiCu with respect to x{9,(f)) and y{9,(()), use {pi,P 2 ,P 3 ) = {.P(,Px,Py) 
and (gi, q 2 , gs) = (C, x, y), and dehne H as 

r ^ ^ 

H = p^uj+ / dE '^piqi-L. (73) 

y i=i 

That procedure yields the same result as the one used above only to hrst order in Gi. 
The terms of second order differ, and the Hamiltonian resulting from Eq. fITHll is not 
rotationally invariant. More precisely: Upon quantization H does not commute with 
the three components Qk of angular momentum in Eq. (|M1) . Therefore, we have not 
used Eq. fl73|) . 


4. Quantized Hamiltonian 

In the present Section we complete the program of the paper. We quantize the effective 
Hamiltonian (ITTjl . we discuss two important discrete symmetries, and we investigate the 
resulting spectra. 


4-1. Quantization 


In quantizing H and the three components Qk of angular momentum, we encounter 
the problem that u and C, are curvilinear coordinates, and that quantization in such 
coordinates is ambiguous, see Ref. and references therein. The quantization depends 
on the physical constraints that limit the dynamics to the curved manifold, and thereby 
on the physical situation. We focus attention on the relevant parts of H and of the Qk- 
These are given by the Hamiltonian of the rigid rotor. 


^rot _ 


2a 




sin^ u 


and by the associated rigid-rotor parts of Eqs. 


(74) 


Q™* = - Puj sin ( - pc^ cot ui cos C , 
Q™* = Plo cos C - Pc cot sin C , 
Qf = PC • 


( 75 ) 
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We use the standard approach [15] to quantizing if™*. In Appendix 4 we describe a 
different approach to quantization which avoids the ambiguity encountered below and 
yields the same result. With 


G-^ = 


1 0 
0 l/sin^o; 

we write ii'’°* in matrix form, 

1 




2C, 




-1 


Puj 

Pc 


(76) 


(77) 


Quantization is achieved upon putting 

1 


^rot _ 


2C'a\/det G 


{—idu}-, —id(f)G ^vdetG 


-iduj 

—idQ 


^d^ + cotu:d^ + 


-cf 

. 2 , 


sm uj 

The transition from Eq. d74|) to Eq. flTHj) is tantamount to putting 

1 


(78) 


Puj= - 


\/sina; 

Pc = - . 


dujVi 


'^Vsincj , 


(79) 


The expressions fl75l) must be symmetrized with respect to prior to using Eqs. fl79|) . 
That gives 

Q™* = ^ sin C i cot u cos ( , 


Q™* = — i cos (dui + 'I cot u sin ( , 


'jrot 
'V 

QT = (80) 

It is easy to check that the three components obey [Q™*, Q™*] = (cyclic). Therefore, 
Q™* = {Q™*, QT^} is ^ bona fide angular-momentum operator. 

For the remaining variables, we impose the usual quantization conditions and choose 
a representation where all the x’s and p’s are ordinary variables so that 


Plu = 


d 


dxi 


pI^ = 


d 


^ l /7 ' dyLfi 

The components of the quantized angular momentum operator are given by 


(81) 


Qx = i sin + i cot u cos C + 


cosC 


sino; 


K , 


sinC 


Qy = — i cos + i cot u sin C "1- ^ > 

sm u 


Qz 1 

with the operator K defined in Eq. 
obey the commutation relations 

[Qx, Qy] = iQx (cyclic) 


(82) 


and the quantization condition flHTl) . The Qk 


(83) 
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The square of the total angular momentum operator is 

— cot H- \ —(—9? + 2ik cosudc + k'^) . (84) 

sin bj ^ 

It is obvious that [Qkik] = 0 for /c = x,y,z. A complete set of commuting angular- 
momentum operators is, thus, Q^,Qz,k. That is expected: For the axially symmetric 
rotor the square of the total angular momentum and its projections Qz onto the 
space-hxed and R onto the symmetry axes are constants of the motion, with quantum 
numbers J(J -|- 1), M, K, respectively. 

We expect that 

[Qk, H] = 0 foT k = x,y,z . (85) 

To prove Eqs. fl85|l we observe that the hrst line of Eq. fITT]) does not involve or 
p(^. Therefore, it suffices to show that the terms in this line all commute with k. 
That is straightforward. We conclude that H consists of two commuting parts: The 
square of the total angular momentum with quantum numbers J{J + 1), M, K, and 
the Hamiltonian for the surface vibrations which carry the quantum number K. As 
a consequence, a rotational band occurs upon every eigenstate of the vibrational part 
of H. All rotational bands have the same moment of inertia. Differences arise only 
through terms not considered in the approximation leading to Eq. fl60|) . see Section ITTI 

4-2. Discrete symmetries 

Discrete symmetries may restrict the spectrum beyond the requirement imposed by 
time-reversal invariance. We follow Bohr and Mottelson [HI [23] and Weinberg [2] . 

We hrst consider 7?.-symmetry. That symmetry is realized if an axially symmetric 
nucleus is, in addition, symmetric under a rotation about vr around an axis perpendicular 
to the symmetry axis. For dehniteness, we choose a rotation r around the y axis. We 
consider the product 5 f( 0 , 6 *)r( 0 , vr, 0 ) where (with operators acting to the right) the 
rotation r is applied prior to g. We have 

5 f( 0 , 6')r(0, TT, 0) = 6, 0)g{0, tt) 

= g{(f) -|- TT, TT — 9)h{7r) . ( 86 ) 

We note that the naive evaluation g{(j), 9)r{0, tt, 0) = gR, 9 + n) would carry 9 outside 
of its domain of dehnition 0 < 9 < n. Eqs. (j 86 il show that the rotational degrees of 
freedom 9 and 0 behave under 7^ as a particle on the sphere under the usual parity, i.e. 
{9, 0) —)■ {n — 9, (f + 7i). The operation TZ acts also on the intrinsic variables x, y, and 
we have 

u{x, y)r{0, TT, 0) = u^—nsintf, Kcos' 0 )r(O, vr, 0) 

= £/(^,«)^(-^)’^( 0 ,vr, 0 ) 

= g{if> + 71,71 — K)h{'if -|- tt) 

= «(—(tt — k) sin('?/^ -|- tt), (tt — n) cos('0 + 7r))h(2'0) 
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= m(—(k — tt) siii'^, (k — tt) cos'ip)h{2f!) 

= u{-x,-y)h{2fj) . (87) 

Here, we used Eqs. fll20p that relate g and u, and geometric considerations in going 
from the second to the third line, and from the sixth to the last line. Eqs. (l87ll show 
that the intrinsic variables transform under TZ as (x, y) {—x, —y)- According to Bohr 
and Mottelson [8l[23], eigenfunctions of the intrinsic and external variables must have 
the same TZ parity, i.e., both must simultaneously be either positive or negative. 

We next consider ordinary parity P. Following Weinberg (Ref. |2], Sect. 19.2), the 
fields X and y have the same parity as the generators Jx and Jy. As components of an 
axial vector, Jx and Jy have positive parity. Thus T'{x,y) = {x,y). That implies that 
all quantized modes (i.e., single excitations, double excitations ...) of the helds x and y 
have positive parity and are allowed. 


4-3. Spectra 


The quantized Hamiltonian is given by Eq. tra with the understanding that the 
momenta are operators as defined in Eqs. d79|) and (IHTj) . The term of leading order 
is given by the first line and written as 



+ (f 


V '|2 


2C, 


+ 


a 2 
2 


2 

Lfi 


+ 2/i/i) 


( 88 ) 


It describes an inhnite set of uncoupled harmonic oscillators with frequencies ojLy. = 
[(L(L + l)Da + In practice, the breakdown scale A serves as a cutoff for 

this sum. It is useful to combine x^y^ and yiy into a two-dimensional S'O (2)-symmetric 
harmonic oscillator with quantum numbers n^y = 0,1, 2 ,..., k^y = 0, ±1, ±2,..., and 
energies {2nLy + {k^yl + The intrinsic angular momentum of the oscillators is 

given by the eigenvalues K = X]l >2 operator K, see Eq. fEB]) . The 

double sum extends over occupied states only. 

For the ground state all quantum numbers vanish. For the excited states, we assume 
Db > 0. The lowest vibrational state corresponds to the single-quantum excitation of the 
mode (x 2 o, l/ 2 o)- As shown in Section IT^ the fields x and y have the same positive parity 
as the corresponding generators Jx and Jy in Eqs. fl36l) . Thus, the lowest vibrational 
state has lA^I = 1 and negative 77.-parity. That is indeed observed in linear molecules |23]. 
Nuclei, however, are different. Here, low-lying states are built from paired Fermions and 
have positive 77.-parity. States with negative 77.-parity correspond to pair breaking and 
have high excitation energies at or above the breakdown scale A of the EFT. Thus, states 
with odd K and positive parity are absent in the low-energy spectrum of nuclei. As an 
example we mention the absence of low-lying magnetic dipole excitations [25l|26], i.e., 
K = 1 states with positive parity, in the spectra of deformed even-even nuclei [271[28| l29]. 

For nuclei, two quanta need to be excited in the lowest (x 2 o, 2 / 2 o) mode, yielding 
a degenerate pair of states with K = Q and [A'l = 2. Data [27l [29] indeed show that 
the low-lying vibrations in even-even deformed nuclei have TT = 0 and lA^I =2. In the 
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actinides and rare-earth nuclei, these states have excitation energies of about 1 MeV, but 
are not degenerate, for two reasons. First, anharmonicities, i.e. higher-order corrections 
to the vibrational Hamiltonian (]88|1 lift any degeneracies, see Section 14.41 Second, the 
proximity of the breakdown scale at A 2 — 3 MeV amplihes such anharmonicities. 
Indeed, there are only few nuclei that exhibit two-phonon excitations on top of the 
one-phonon K = 0 ot \K\ = 2 vibrations, see Refs. [301 El] for recent reviews. The 
impressive spectra of ^^^Er [27] and of ^®^Dy [29] conhrm this picture. The positive- 
parity states in those spectra must be viewed as anharmonically distorted quantized 
vibrations corresponding to our Nambu-Goldstone modes. These spectra also show 
negative-parity states. These states cannot be understood within the EFT discussed in 
this paper; they can possibly be viewed as vibrations on top of the intrinsic odd-parity 
state with lowest energy. 

The quantized version of the full Hamiltonian of Eq. (ITT!) is 

J 2_^2 




2Ca 


(89) 


The last term causes a rotational band to appear on top of each of the vibrational 
states (band heads). The eigenfunctions are Wigner D-functions j^{a, (3,0) with 
total integer spin I and projections —I < M, K < I usms]. The eigenvalues of P are 
/(/ -|- 1) with I > \K\. At this order, all rotational bands have the same moment of 
inertia, and deviations from this picture are due to higher-order corrections, see Ref. [TO] 
and Section 031 

We conclude that the EFT predicts that in leading order, the Nambu-Goldstone 
modes due to emergent breaking of rotational symmetry yield a large number of 
harmonic vibrations. In next-to-leading order each of these serves as head of a rotational 
band. All bands have identical moments of inertia. Gorrections of higher order 
considered in Section 14.41 lead to anharmonicities of the vibrational states and cause 
the moments of inertia to differ. 

Given the close proximity of the breakdown scale in nuclei to the vibrational 
excitation energy, it is reasonable to consider a simpler - but still model-independent - 
approach to deformed nuclei. That approach [15] uses an effective theory (as opposed 
to the effective held theory of the present paper). It combines the quantized rotations 
as lowest-energy excitations with the lowest vibrational modes. The latter correspond 
to the K = 0 and K = 2 modes of the present paper. Thus, the effective theory replaces 
the quantum helds x and y by their quantized modes of longest wave length. Our results 
show that the effective theory is based upon a solid held-theoretical foundation. 


f.f- Terms of Higher Order 

The extension of the effective held theory to higher-order terms is straightforward 
but tedious. The program is this, (i) Use the power counting of Section 13.41 to 
identify all terms that contribute to the ehective Lagrangian in a given order. That 
includes terms with higher time derivatives. These can be treated by perturbative 
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field redefinitions [321 EH]- (ii) Expand the fields x and y into their normal modes as 
in Eqs. fl59|) . and compute the Lagrangian by integration of the Lagrangian density 
as in Eqs. (|57|l . The resulting expressions are (complicated) sums involving Clebsch- 
Gordan coefficients, (hi) Perform the Legendre transformation to the Hamiltonian 
within perturbation theory to the desired order of the power counting, (iv) Quantize 
the Hamiltonian and compute the spectrum. 

Obviously, steps (i) to (hi) are quite laborious. Furthermore, the kinetic part of 
the resulting Hamiltonian will have low-energy constants that are complicated linear 
combinations of the corresponding coefficients of the Lagrangian. The latter are not 
known, and the former need to be determined from data. It is, therefore, desirable 
to understand the transformation properties of the coordinates xl^ and yi^ and of the 
canonical momenta and and to directly construct the most general Hamiltonian 
that is invariant under rotations at a given order of the power counting. 

The construction of the invariants is guided by the following observations. Eqs. 6D 
show that under rotations, the helds x{9, cf) and y{9, (f) transform as the x and y 
components of a two-dimensional vector. These transformation properties hold for 
every point (0, 0) on the unit sphere. Using the expansion of the second of Eqs. (15^ 
we conclude that the complex normal modes xl^ and yi^ themselves, too, transform 
as the X and y components of a two-dimensional vector. The same is true for the 
corresponding canonical momenta (denoted by pf^ and respectively) as these stem 
from time derivatives of the held modes xl^ and y^^. The integration J dE over the 
Lagrangian density in Eq. (1571) singles out scalars. For instance, the invariant that is 
bilinear in the momenta is 


-pin Pi- Pi) = fcPL-M + PlpPl-^) 

L Lfi 


= E {(piA + (plA) ■ 

L{i 


(90) 


Here pf and pl denote spherical tensors of degree L with components pf^ and pl^, 
respectively. 

We apply these considerations hrst to the kinetic terms in Eqs. (154)1 and (155|1 
and then to the potential terms in Eqs. fl56|) . In calculating the kinetic part of the 
Hamiltonian we encounter the need to invert the generalized form of the matrix G in 
Eq. fl70|) . So far we have taken into account only the terms Cia and Cib in Eqs. (154)1 . 
As in Eqs. ()57)l we now consider the sum £i of all four terms in Eqs. ()54)l and ()55)l . 
The kinetic part Li of the effective Lagrangian is given by integration over dE of Ci. 
We note that Ci is bilinear in the time derivatives of all dynamical variables. We omit 
the term which gives a trivial contribution. Dehning Xj as the totality of the time 
derivatives we proceed as in Section 1331 and write 


Li = 



(91) 


In order not to overburden the notation we have chosen the letter A rather than G for 
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the matrix defining Li. As in Section 13.51 we write A for the matrix and Aij for its 
elements. We have Aij = Aji, and we write A = + A^^\ Here is the 

sum of the contributions arising from Cia and Cih (from Cic and Cid, respectively). The 
latter are small of order £ with respect to the former. The momenta are defined by 
Pi = AijXj, and the kinetic part of the effective Hamiltonian is 

^1 = • (92) 

ij 

We calculate A~^ perturbatively and use the explicit form of Cic and Cid in Eqs. fl5Tl) 
and ([55]). Details are given in Appendix 5. The contributions of order De to the effective 
Hamiltonian stemming from Cic (-^id) are denoted by Hic {Hid, respectively). We find 

a 


Hic= - 


Hid= - 


2Cl 

2Cl 


dD {xpy - ypA^ , 

dD {x"^ + yAipl + pI) 


A contribution of order De arises also from Cu- It is given by 

1 COS^ OJ 


Hh = 


dfi {xpy - ypxY 


(93) 


(94) 


2Ca sin^ uj 

We turn to the potential terms. For the modes xl^. and yLy, we deal, in analogy to 
Eqs. ([Mil , with the invariant ~kyLii)- evaluating the remaining invariants 

in Eqs. fl56|) we have to deal with the angular-momentum operators L and Lz acting on 
the fields x{9, 0) and y{9, 0). Denoting by Lip the spherical components of L, we have 


(95) 


LyX{0,4’) ^ \LyYly{0,4>) 

Lfi 

Lfi 

Upon multiplication with Yf^{9,(j)) (where a and a are arbitrary) and integration over 
dE we find 

Lfi 


AEY^»,4,)Px(e,i,) = 


AEY:„(e,4')YLy+y(»,4') 


— no-ct. 


'luLa-v^aa-u 

= {-lYL_,pXaa ■ (96) 

In the last line the operator L_ip is understood to act on the spherical tensor of rank a 
with components Xaa- Invariants built upon the normal modes xl/i and ym are, thus, 
obtained by viewing these modes as components of spherical tensors of rank L. Scalars 
are constructed after acting with L onto these components. For example 


jAE (Lxf = ^(-1)" j 






aa by 

aa by 


V 

a+ix 
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aa 


^aa^a—a 


= - XL- XL • 


(97) 


The evaluation of the terms in Eqs. fl56l) is now straightforward. Upon applying the 
quantization rules, we obtain higher-order terms in the Hamiltonian H. 

The terms in Eqs. (jMl) and (j56il do not depend on the rotational degrees of freedom 
u and (. These terms lift the degeneracies of the vibrational modes of excitation but 
do not affect the moment of inertia. The term in Eq. fIMll depends parametrically upon 
u and couples the rotational bands with the vibrational modes. 


5. Summary and Conclusions 

We have constructed the EFT for emergent symmetry breaking in deformed nuclei. In 
addition to the Nambu-Goldstone modes, the theory contains two additional degrees 
of freedom that describe rotations about the body-fixed symmetry axis. Starting from 
a physically intuitive and mathematically standard parameterization where rotational 
and vibrational degrees of freedom are treated on an equal footing, we have switched to 
a much more practical parameterization where the rotational degrees of freedom receive 
a separate treatment. The theory is characterized by three small parameters. These 
are (i) the ratio of the energies of rotational motion and of vibrational motion 
H, (ii) the ratio H/A where A is the breakdown scale of the EFT (typically given by 
the pairing energy or the energy of single-particle excitation in the shell model), and 
(hi) the parameter e which characterizes deviations from harmonicity of the vibrations. 
In lowest order, the spectrum consists of vibrations each of which serves as head of a 
rotational band. The vibrations are due to the quantized Nambu-Goldstone modes that 
describe the emergent breaking of SO(3) symmetry. In leading order, the vibrational 
modes are degenerate, and the rotational bands all have the same moment of inertia. 
Terms of next order remove both degeneracies. 
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Appendix 1: Equations of Motion in 
Curvilinear Coordinates 

In curvilinear coordinates the equations of motion are obtained by variation of 
the product of the Lagrangian density C and the integration measure. We 
demonstrate that fact for the simplest case. We consider a Lagrangian density 
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Cigp^dfj/dt^dfj/dxi,... /dxjsi) that depends on the held (a function of time t 
and of N Cartesian variables Xi,X 2 , ■ ■ ■, xjsf), and on the iV + 1 derivatives of the held. 
Standard variation of the action integral yields 

^ ^ N O N 


dt 


JJdx, 


d dC 




d dC dC 


+ 


U =1 


dtdidtfj) ^ dXf,d{d^ fj) difJ 

LL —i 


■ 


(98) 


We introduce N curvilinear coordinates (k{xi,X 2 , ■ ■ ■ ,xn) with k = 1,..., iV that are 

irte 

Then 


functions of the N Cartesian coordinates We dehne the iV-dimensional matrix 
M,,k = with D = det M . 


dCk 


(99) 


N 


N 


dx,, = D dCfc , 

u=l k=l 

We insert all this into expression 


d 


N 


dx 




d 


dC 


N 


^ k=l 

and obtain 


dCk ’ 5(4 ) 


rtr 

-( 100 ) 


1=1 


5(4) 


N 


dt 




5 dC 


k=l 

N N f) ^ 

/i=i 1=1 


dt d{dt'ip) 


n=l 


dC dC 


)H’ 


( 101 ) 


The triple sum in expression (llOlh can be written as 


N 

E 


5 dC 

^ dCi 5(5^4) 




( 102 ) 


1=1 \ sf / filn=l 

The identities det In D = In Trace D and = d^^M^i imply that expression (11021) 

is equal to 


5 dC 


N 

^ dCi 5 ( 5 ^ 4 ) 


D 


-1 


dC 


y—D, 

d(, I d(d^^i<) 


(103) 


Using all that and the independence of H of t and ip we rewrite expression fllOip as 


„ N 

dt / JJdCfc 

k=i 


+ E 


5 d{DC) d{DC) 


5 d{DC) 

dt didtip) ' ^ dCi didQ'ip) ' dip 


+ 


5 ip 


(104) 


Comparing expression fll04p with expression fl95p we conclude that variation of C 
in curvilinear coordinates is tantamount to varying DC and otherwise treating the 
curvilinear coordinates like Cartesian ones. That is what we use in Section 12.31 with 
D = sin6'. 
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Appendix 2: The Matrix M 


We calculate the matrix M defined by 5qi, = 'YlkMuk^Xk- We define 
H = u [cos (Ex' + sin (Py'] 5 

G = oj [—sin(^Pa./+ cos^Py/] . (105) 

Then 

[P, G] = , [P, -iP,,] = G . (106) 

With r given by Eq. (|23|1 we write 

rexp{-iP} = exp{-zP}exp{-i5^ P^,} . (107) 

Here 

H = Cj [cos C Px' + sin C, Py'] , 

OJ = u + 5oj , 

C=C + 5C. (108) 

We calculate Su,S( and 6 ^ to hrst order in Sxk- Keeping only terms of hrst order we 
have with k = x',y',z' 

dxkPk = —P + Sf, exp{—iH}Pz' exp{iH} + 6(X (109) 

k 

where 

00 / .\^_2 k —1 

A'= (X^rZ^j—5^tfG//'=-‘-‘)exp{i//} . (110) 

fc=0 ■ 1=0 

Since X must be of order zero in P only terms with / = fc — 1 in Eq. (I110|l contribute. 
From the commutation relations fll06|) we have to zeroth order in P 



H^G —)■ Gu’’ for 1 even , 

H^G - iP,,u^+^ for 1 odd . 

(Ill) 

Thus, 

X = (— sin^Pj,/ + cosCPyO sino; + P 2 /(cosa; — 1) . 

(112) 

Similarly, 

H\—iPzi) —)■ Gu^~^ for 1 odd , 

H\—iPz>) — iPz'OJ^ for 1 even . 

(113) 

Therefore, 

exp{—iPjPj,/exp{iP} = P^' cosa; + {—sm(P^i + cosCPy) sinca . 

(114) 

Insertina; Eqs. flll4p and flll2p into Eq. fllOOp we hnd 



SxkPk = du}{cos CPx' + sin (Py,) - P^^SC 



+ + SC)[Pz' coscn + {—sin^P^i + cos^P^/) sina;] . 

(115) 
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We equate the coefficients multiplying Pyi, and P^/ on both sides of that equation. 
With 6 u = 6 qi, 6 ( = 6 q 2 , 5.^ + (5C = we obtain 

3 

SXk = ^{M-^)ki5qi (116) 

1=1 


where 


M-^ = 

The inverse matrix is 


/ cosC 
sinC 

V 0 


0 

0 

-1 


M 


cos C 

sin (f cot oj 

sin C 


— sin C, sin u' 
cos C sin u 
cos UJ 


sinC 

cos C cot OJ 

cos f 
sincj 



(117) 


(118) 


Appendix 3: Equivalence of the 
Parameterizations (I36p and ([4], [5]) for U 

We start from Eqs. fl36l) and derive Eqs. (jU [5]). We dehne k > 0 and T by 
X = — K, sin T , 

y = Kcos T . (119) 

We use a power-series expansion in k, valid for k -C 1. We prove the equivalence only 
to leading order in k. Terms of higher order can be treated analogously. In addition to 
Eqs. fl36|l we use the following dehnitions and identities, valid for arbitrary values of k, 
of the Euler angles a, (3,'y, and of T, 

h(7) = , 

r{a,/3,'y) = g{a,/3)h{-f) , 

m(—K sinT, KCOST) = 5 f(T, K)h^(T) . ( 120 ) 

The hrst of Eqs. (I12np dehnes h{'j) for an arbitrary angle 7 . The second of Eqs. (11201) 
dehnes the rotation r{a,(3,'j) and shows that in Eqs. fl36|l . the factor g{C)^) ^-cts on u 
like a rotation r with third Euler angle 7 = 0 . The third of Eqs. (11201) is an identity 
for the function u{x,y) dehned in the third of Eqs. (|36ll . That identity can easily be 
derived with the help of the relation exp{—iT7,^} exp{iT7,^} = — sinTJ^, -|- cosTJ^. 

We use the transformation law 

^(C, 0 )^(T, k) = ^(C', uj')h{y) . ( 121 ) 

Here (',oj','y' are functions of the angles ( and u and of the variables T and k and are 
given in Ref. [T5|. We have 

sincj 

cot(C — 0 = cos a; cot w -|- cot k ——, 

sm W 

cos oj' = cos K cos (u — sin K sin oj cos T , 

, sin K 

cot 7 = — cos K cot W — cot a;-— . 

smT 


(122) 
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For |k| -C 1 that yields in leading order 


C' 

_y 


~ a; + K cos ~ o; + 1/ , 

^ Ksind' ^ X 

~ c + —— ~ c - ^—, 

sin uj sin uj 

~ 'F — K cot a; sin d' = d' + a; cot u . 


(123) 


Therefore, 


^ g{C - + y)h{x cot u) . 

smuj 


(124) 


We use the third of Eqs. 


once again to write the result as 


/ ^ X \ 

^u[- {u + y)sm{C-- -), (a; + i/) cos(C - ^-) 

V sincj sincj / 

X h{'^)h{xcotu) . (125) 


We recall that U is dehned in the coset space S03/S02. Hence 

U ^u(-{u: + y) sin(C - .{u + y) cos(C - ^^)) • (126) 

V sincj sincj / 

Eq. (I126p gives the connection with the parameterization of U in Eqs (01) and (j5]) to 

lowest order in x, y. Differences in sign are due to the fact that here we work in the 

space-hxed system. The variables ojq and Co in Eqs. (jS]) are seen to correspond to u 

and Cl respectively. To lowest order in k, the variable cji corresponds to y whereas 

the variable Ci corresponds to —x/ sino;. The occurrence of the factor 1/ sinca in the 

last relation causes the difficulties in the attempt to derive the equations of harmonic 

motion directly from the parameterization (0]) via an expansion in powers of uji and Ci 

and explains why we have introduced the parametrization fl36l) . The calculation can 

obviously be carried to higher orders. That establishes the complete equivalence of the 

parameterizations of U in Eqs. (|36ll and (jUlS]). 


Appendix 4: Quantization in Curvilinear 
Coordinates 

It is instructive to use another approach to quantization which shows how the 
ambiguities that are associated with the prescription fl79|) are avoided. We assume 
that oj and C are independent so that 

[n;,C] = 0 , = 0 , [p^,C] = 0 , [p^,Pc] = 0 • (127) 

For simplicity we have suppressed the symbol " on the operators. It remains to 
determine the commutators [pu},oj] and [pc,C]- We choose a representation where oj 
and ( are ordinary real variables. Quantization is subject to three requirements, (i) 
The expressions for Qy* and for if™* must be Hermitian. (ii) The Qy* must obey the 
standard commutation relations (cyclic), (iii) When expressed in 

terms of the quantized components Q™* of angular momentum, the Hamiltonian of the 
pure rotor must be given by Eq. (I35D . with C —)• Ca- As for point (i), Hermitecity is 
dehned with respect to an integration measure for the variables oj and Q The matrix U 
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introduced in Eqs. (I5B]) is defined in the coset space SO{3)/SO{2). That suggests using 
for the volume element the expression 


(128) 


dV = do; sin a; dC • 

An operator O is Hermitian if the equality 



(129) 


holds for any two integrable functions ^(a;, () and T(a;, () that are periodic with period 
27r with respect to both uj and (. (That is why we have chosen in Section [3] the ranges 
of integration as 0 < uo, Co < 27r). Requirement (ii) on the commutators of the Q™* and 
the explicit form of the Q™* in Eqs. fl7^ imply 


= -i , |p<.cl = -i ■ 


(130) 


Eqs. (1751) show that for to be Hermitian, and p(^ must be Hermitian, too. As 
formulated in Eq. fll29p . that condition is consistent with Eqs. fll30p if p^ and obey 
Eqs. fl79l) . Calculating from Eq. fl35ll with C —)■ Ca we obtain the second of 
Eqs. fITSp . Thus, if™* is the quantized Hamiltonian of a rotor with Ca the moment of 
inertia. 

Appendix 5: Kinetic Part of the Effective Hamiltonian 


We use i-* ^ (i(o))-* - (i(o))-*iW(i(o))-* + .... With 



(131) 


we have 



(132) 


The non-diagonal terms cx are of order Therefore, we invert by 

expanding in powers of these terms and obtain to first order 



(133) 


The matrix receives contributions from both Cic and Cid- The contribution from 
Lie is proportional to 



LL'L"L"'Li 



( 134 ) 













Effective Field Theory of Emergent Symmetry Breaking in Deformed Atomic Nuclei 32 


The upper indices on the big round and square brackets denote the total angular 
momentum to which the terms are coupled. For Li^ we have correspondingly 




in 0 


LL'L"L’"Li 

+ 2 Q cos oj 


Y. 

LVL"L"'Li 


XLfiXL'fi' + yL^lyL>^l> 


Li 


X I XL"^_L'|yL"'^^"| — XL"|J,"yL"'^l"' 


Lr 


(135) 


Both Lie and Li^ are small. Therefore, we calculate —to lowest 
order, i.e., by taking for only the diagonal part of the matrix on the right-hand 

side of Eq. (11331) . The resulting terms in the effective Hamiltonian are then given by 


Hu = 


Hm = 


a 


E 


2C? 

^ LL'L''L"'Li 


Li 


XlfiPiifj^/ PLij.yL'iJ.' 

, in 0 

X \XL"ii"Pf infill - PL"ti"yL"'/i"'' 


Ca 

2Cl 


E 


ii'i''i"'ii 


XL^lXL'^^' + yL^iPL'fj.' 

Lin 0 


X \PL"^"PL"'^l"' + p\"^_ l'iP\iii fj,'" 

Cd 2 pQ COS OJ 


2 CaCb sin^o; 


E 


LL'L"L’"Li 


XLfiXL'ii' + yLtiPUti' 

Li"! 0 


Li 


X L"bt''P\i"fj,'" PL"bi"yL"'fi"''^ 


(136) 


The terms inversely proportional to are of order eQ, the term inversely proportional 
to CaCb is of order and, therefore, negligible. That result justifies a posteriori the 
diagonal approximation for With the help of the expansions fl5^ and fl6^ . we 

obtain Eqs. fl^ . 
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